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Abstract 

In experiments with the CMD-2 detector at the VEPP-2M electron-positron 

collider at Novosibirsk about 150000 e^e~ tv'^tv^ events were recorded in the 
center-of-mass energy range from 0.61 up to 0.96 GeV. The result of the pion form 
factor measurement with a 1.4% systematic error is presented. The following values 
of the p-meson and p — u; interference parameters were found: Mp — (775.28±0.61± 
0.20) MeV, Tp = (147.70 ± 1.29 ± 0.40) MeV, r(p ^ e+e^) = (6.93 ± 0.11 ± 0.10) 
keV, Br{u! ^ tt+tt") = (1.32 ± 0.23)%. 
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1 Introduction 



The cross-section of the process e+e tt+tt is given by 

a^—Pl\FAs)\\ 

where -FV(s) is the pion form factor at the center-of-mass energy squared s 
and /Jjr is the pion velocity. 

The pion form factor measurement is important for a number of physics 
problems. Detailed experimental data in the time-like region allows mea- 
surement of the parameters of the p(770) meson and its radial excitations. 
Extrapolation of the energy dependence of the pion form factor to the point 
s = gives the value of the pion electromagnetic radius. Exact data on the 
pion form factor is necessary for precise determination of the ratio 

R — a{e^ — > hadrons) / a{e^ ^'^ jjT). 

Knowledge of R with high accuracy is required to evaluate the hadronic con- 
tribution a^'^'' to the anomalous magnetic moment of the muon [g — 2)^ 
About 87% of the hadronic contribution in this case comes from s < 
2 GeV^/c^ (VEPP-2M range), and about 72% — from e+e" tt+tt" chan- 
nel with s < 2 GeV^/c^ g ^. The E821 experiment at BNL g has collected 
its first data in 1997 and will ultimately measure {g — 2)^ with a 0.35 ppm ac- 
curacy. To calculate the hadronic contribution with the desired precision the 
systematic error in R should be below 0.5%. Therefore a new measurement 
of the pion form factor with a low systematic error is required. 

Experiments at the VEPP-2M collider[||, which started in the early 70s, 
yielded a number of important results in e~^e~ physics at low center-of-mass 
energies from 360 to 1400 MeV. The high precision measurement of the pion 
form factor at VEPP-2M was done in the late 70s - early 80s by OLYA 
and CMD groups |§. In the CMD experiment, 24 points from 360 to 820 
MeV were studied with a systematic uncertainty of about 2%. In the OLYA 
experiment, the energy range from 640 to 1400 MeV was scanned with small 
energy steps and the systematic uncertainty varied from 4% at the p-meson 
peak to 15% at 1400 MeV. 

During 1988-92 a new booster was installed to allow higher positron cur- 
rents and injection of the electron and positron beams directly at the desired 
energy. During 1991-92 a new detector CMD-2 was installed at VEPP-2M, 
and in 1992 it started data taking. The pion form factor measurement was 
one of the major experiments planned at CMD-2. The energy scan of the 
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ID 


Date 


2E, GeV 


Number of 
energy points 




1 


Jan-Feb 1994 


0.81-1.02 


14 


35000 


2 


Nov-Dec 1994 


0.78-0.81 


10 


66000 


3 


Mar- Jim 1995 


0.61-0.79 


20 


85000 


4 


Oct~Nov 1996 


0.37-0.52 


10 


4500 


5 


Feb-Jun 1997 


0.98-1.38 


37 


75000 


6 


Mar-Jun 1998 


0.36-0.97 


37 


1900000 



Table 1: CMD-2 runs dedicated to R measurement 



whole VEPP-2M 0.36-1.38 GeV energy range was performed in six separate 
runs listed in Table |l|. The energy range below the 4> nieson was scanned 
twice (in 94-96 and in 98). 

In this article we present results of the analysis of the data from runs 
1-3. The data were taken at 43 energy points with the center-of-mass energy 
from 0.61 GeV up to 0.96 GeV with a 0.01 GeV energy step. The small 
energy step allows calculation of hadronic contributions in model-independent 
way. In the narrow energy region near the w-meson the energy steps were 
0.002 0.006 GeV in order to study the w-meson parameters and the p — ui 
interference. Since the form factor is changing relatively fast in this energy 
region, it was important that the beam energy was measured with the help 
of the resonance depolarization technique at almost all energy points. That 
allowed a significant decrease of the systematic error coming from the energy 
uncertainty. 

The CMD-2 (Fig. % is a general purpose detector consisting of the drift 
chamber, the proportional Z-chamber, the barrel (Csl) and the endcap (BGO) 
electromagnetic calorimeters and the muon range system. The drift chamber, 
Z-chamber and the endcap calorimeters are installed inside a thin supercon- 
ducting solenoid with a field of 10 kGs. More details on the detector can 
be found elsewhere |0, ^, ||]. The data described here was taken before the 
endcap calorimeter was installed. 

Two independent triggers were used during data taking. The first one, 
charged trigger, analyses information from the drift chamber and the Z- 
chamber and triggers the detector if at least one track was found. For 
0.81-0.96 GeV energy points there was an additional requirement for the 
total energy deposition in the calorimeter to be greater than a 20-30 MeV 
threshold. The second one, neutral trigger, triggers the detector according 
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Figure 1: Detector CMD-2. 1 beam pipe, 2 drift chamber, 3 — Z- 
chamber, 4 — main solenoid, 5 — compensating solenoid, 6 — endcap (BGO) 
calorimeter, 7 — barrel (Csl) calorimeter, 8 — muon range system, 9 — 
magnet yoke, 10 — storage ring lenses 

to information from the calorimeter only. Events triggered by the charged 
trigger were used for analysis, while events triggered by the neutral trigger 
were used for trigger efBciency monitoring. 

2 Data analysis 

2.1 Selection of collineair events 

The data was collected at 43 energy points. At the beam energy of 405 

MeV the data was collected in two different runs with different detector and 
trigger conditions. Therefore, two independent results for this energy point 
are presented. 

From more than 4 • 10^ triggers about 4 • 10^ were selected as coUinear 
events. The selection criteria were as follows. 

1. The event was triggered by the charged trigger. Other triggers may be 
present as well. 
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Figure 2: Definition of parameters for two track events 

2. Only one vertex with two oppositely charged tracks was found in the 
drift chamber. 

3. Distance from the vertex to the beam axis p is less than 0.3 cm. 

4. Z-coordinate of the vertex (distance to the interaction point along the 

beam axis) \Z\ is less than 8 cm. 

5. Average momentum of two tracks {p\ -\-p2)/2 is between 200 and 600 



6. AcoUinearity of two tracks in the plane transverse to the beam axis 
|A(^| = Itt — |yi — <^2|| is less than 0.15 radians. 

7. AcoUinearity of two tracks in the plane that contains the beam axis 
I AO I = |0i — (tt — ©2)1 is less than 0.25 radians. 

8. Average polar angle of two tracks [Oi + (tt — 02)]/2 is between 0mm 
and (tt — Qmin)- AH analysis was done separately for Qmin = 1-0 and 
©mm = 1-1 radian. 

Because of the unfortunate accident the high voltage was off for two neigh- 
bouring lines of the Csl calorimeter during the data taking in the 0.61 0.784 
GeV energy region. In order to avoid additional calorimeter inefficiency, 
events with 4.35 < (^+ < 4.95 or 3.90 < < 4.50 have been removed 
from analysis for corresponding energy points (angles are measured in radi- 
ans). Events with 1.05 < < 1.65 or 0.55 < < 1.15 were additionally 
removed for 0.65-0.70 GeV energy points. 



MeV/c. 
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Definition of the coUinear event parameters, such as p, 61,2, 'Pi,2j is iUus- 
trated in Fig. ||. 
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2.2 Event separation 
2.2.1 Likelihood function 

The selected sample of events con- 
sists of collinear events e'^e^ e^e^, 
e~^e^ — > TT+TT", e'^e^ — > fi'^fJ.^ and 
background events mainly due to cos- 
mic muons. The energy deposition in 
the barrel Csl calorimeter by negati- 
vely {E~) and positively {E~^) charged 
particles were used for event separa- 
tion. The distribution of E'^ versus 
E~ is shown in Fig. ^. Electrons and 
positrons usually have large energy de- 
position since they produce electromag- 
netic showers. Muons usually have small 
energy deposition since they are mini- 
mum ionising particles. Pions can in- 
teract as minimum ionising particles 
producing small energy deposition or 
have nuclear interactions inside the calo- 
rimeter, resulting in long tails to a higher 
value of the energy deposition. 

The separation was based on the minimization of the following unbinned 
likelihood function: 



200 



400 



E-, MeV 



Figure 3: E^ versus E~ distribu- 
tion for collinear events for 0.8 GeV 
energy point 



(1) 



where a is the event type (a — ee, /i/i, tttt, cosmic), Na is the number of 
events of the type a and fa{E'^, E^) is the probability density for a type a 
event to have energy depositions E^ and E^ . It was assumed that E^ is 
uncorrelated with E~ for events of the same type, so we can factorize the 
probability density: 



UE+,E-) = f:{E+)-f-{E-). 
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For e^e^, fJ-^ pairs and cosmic events the energy deposition distribution 
is the same for negatively and positively charged particles, while these dis- 
tributions are significantly different for tt"*" and 7r~. Therefore = f~ for 
a = ee, ii^ and cosmic, but for pions these functions are different. The 
probability density functions are described in detail in the following sections. 

In minimization the ratio -/V^^/A'^ee was fixed according to the QED cal- 
culation 

Nf,f, _ Q-pp ■ (1 + Sf^^) (1 + g^^) 

Nee fee ' (1 + ^ee) (1 + 

where a is the Born cross-section, (5 is a radiative correction, a is the cor- 
rection for the experimental resol utio n of O angle measurement and e is the 



reconstruction efficiency (see sec. 2^ for details). The likelihood function (|l|) 



was rewritten to have the following global fit parameters: 

{Nee N^fj^), ~jr^ T,j I Neosmic 

instead of Nee, N^^^^, TV^^r and N cosmic- The number of cosmic events Neosmic 
was determined before the fit as described below, and was fixed during the 
fit. 



2.2.2 Rejection of background events 
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Figure 4: To determination of num- 
ber of cosmic background events 



Cosmic background events can be 
separated from the beam-produced col- 
linear events by the distance p from 
the vertex to beam axis and (or) by 
the Z-coordinate of the vertex. The p- 
distribution for coUinear events has a 
narrow peak around 0, while the same 
distribution for the background events 
is almost flat (Fig. 5(a)). Similarly, 
the Z-distribution for coUinear events 
is Gaussian like, while it is almost fiat 
for the background events (Fig. |5(b)| ). 
Assuming that p- and ^'-distributions 
for cosmic events are not correlated 
and that all events with p > 0.5 cm 
or \Z\ > 8 cm are cosmic events, one 
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Figure 5: p- and Z-distributions for background and coUinear events. Empty 
histograms correspond to aU selected events, filled histograms correspond to 
background events only 



can calculate the number of background events in the signal region from the 
number of events out of the signal region. 

Let us define Nx as the number of cosmic events in the signal region 
(p < 0.3 cm, \Z\ < 8 cm) and iVi.2,3 as the number of events in different 
out-of-signal p — Z regions (see Fig. 

p<0.3; 8<|Z|<12 
0.5<p<1.0; \Z\<8 
0.5<p<1.0; 8<|Z|<12 

Then under assumptions mentioned above we can derive: 

Nx ^ — ■ Ni = — ■ N2. 




Nx is equal to Ncosmic only statistically. The confidence interval for TV, 
corresponding to one standard deviation, is 



COS7mC7 



Nx±V5N^. 
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In order to take that into account, the N cosmic value was fixed to the Nx 
value during fit, but AA'x was added to AA'^jr at the end. 



2.2.3 Energy deposition parametrization 

In order to construct the likelihood function (|^) it is very convenient to have 
an analytical form for energy deposition distributions. We parametrize all 
energy depositions by the linear combination of two types of Gaussian-like 
functions, described below. 



Normal distribution g{x\XQ^a) 
g{x;xo,cr) = 



27rcr 



exp [-{x-xof/2a^] , 



where the parameters are the mean value xq and the standard devi- 
ation (7. The Gaussian is normalized at the (— cx),+oo) interval, but 
energy deposition values belong to the [0, -l-oo) range. Therefore for 
parametrization the Gaussian normalized at the [0, -l-oo) range is used: 



gix;xo,a)/ 



where 



I{x) 



1^ 



exp 



dt. 



Normal logarithmic distribution[10| gi{x;xQ,a,r]) 



gi{x;xo,(T, rf) 
where 



1 



1 



2iTa ctq 



• exp • 



In^ (1 - 2--^7?) 



fJo = In (?7V21n2 + ^/lT^^f'^2]^ j %/2 In 2 . 

Parameters of this function are the most probable energy xq , the a = 
FWHM/2.35 and the asymmetry r\. If 77 = 0, the logarithmic Gaussian 
is equal to the plain Gaussian: 

gi{x\ xqg, f]) \r,=o = g{x; xq, a). 

The function is normalized at the [— 00, xo -|- cr/77] interval: 

rxo+a/ri 

gi{x)dx = 1, 
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but in our case its integral over the range (~cxd, 0] is negligible and it is 
assumed that this function is normalized at the [0,xo + a/rj] interval. 

Since both functions are normalized, it is easy to construct the normalized 
energy deposition distribution. 



2.2.4 Correction of energy deposition for polar angle 

Energy deposition in the calorimeter is correlated with the polar angle of a 
track due to the dependence of calorimeter thickness on 8. So the probability 
density for an event of type a to have energy deposition E depends not 
only on E, but also on Q. In order to simplify the likelihood function, the 
energy deposition was corrected to make dependence of energy deposition on 
Q negligible. In other words, instead of fa{E, 9) we are using fa{E), where 
E is the corrected energy deposition: 

E^E to have fa{E,Q) ^ fa{E). 

The correction procedure was performed as follows. First, two groups 
of events were selected: cosmic events and electrons. Selection criteria for 
both types of events are described below in sections ^.2.5 and 2.2.6, The 



scatter-plot of the energy deposition versus O for coUinear events is shown 
in Fig. |6^. The full range of possible values [1, tt — 1] was divided into 50 
equal bins and for each bin an average energy deposition was calculated for 
electrons and cosmic events. 

The energy deposition of electrons depends on 6 because the CMD-2 
calorimeter has only about 8Xo and therefore the shower leakage is not neg- 
ligible. The effective calorimeter thickness increases for incident track angle 
away from normal. The corresponding small variation of the average energy 
deposition of electrons was fit to the parabolic function 

Ee{Q)^Ei- [1 + a(e - 7r/2)2] , 

where Ei and a are the fit parameters. 

Cosmic events contain minimum ionising particles and their energy de- 
position is proportional to the length in the calorimeter through which the 
particle passes. This length is proportional to l/sin0 , and the fit function 
in this case was 

E,{Q)^Eo/BinQ, 

where Eq is a fit parameter. An example of the fit of the average energy 
depositions for electrons and cosmic tracks is shown in Fig. 
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(a) Energy deposition versus 
for collinear events before cor- 
rection 
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(b) Average energy deposition 
versus © for electrons and cos- 
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correction 
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(c) Energy deposition versus 
for collineax events after cor- 
rection 



(d) Average energy deposition 
versus for electrons and cos- 
mic background events after 
correction 



Figure 6: Correction of energy deposition for the polar angle of the track 
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The corrected energy was calculated as 



E[E,Q) = E- 



where fcc(e) = E^/E^{Q), fce(e) = Ei/Ee{Q), EdO) and E^iQ) arc the 
measured average energy deposition for cosmic and Bhabha events respec- 
tively. After such a correction the average energy deposition for both cosmic 
and Bhabha events does not depend on 8. In fact, ii E ^ Ee{Q), then 
E = Ei; a E ^ Ec{Q), then E = Eq. 

Since the energy deposition of cosmic events does not depend on the beam 
energy, the variations of Eq reflects residual variations in the calibration of 
the Csl calorimeter. To correct for these variations, the energy deposition 
was normalized to have the average energy deposition of the cosmic events 
equal to 85 MeV (the arbitrary constant close to the experimental value) for 
all beam energies: 

- - 85 
E^E-—. 

Eo 

The scatter-plot of the corrected energy deposition versus Q for coUinear 
events is shown in Fig. The average energy deposition for electrons and 
cosmic events versus Q after the correction is shown in Fig. ^d. 

The described energy deposition correction was calculated independently 
for each energy point. In the all following sections E actually means the 
corrected energy E. 



2.2.5 Energy deposition of cosmic events 

As mentioned above, it is possible to select a pure sample of background 
events from the real data. In this case the selection criteria are the same as 
standard cuts for the coUinear events, except that the distance p from the 
vertex to the beam axis should be from 5 to 10 mm (instead of less than 
3 mm). In order to make the selected sample of background events even 
cleaner, the fact that there are two clusters in the event was used. Analysing 
the energy deposition of one cluster, the energy deposition of the other one 
was required to be less than 150 MeV. Since the clusters are independent, 
the strict limit on one cluster energy does not affect the energy deposition 
distribution for the other cluster. 

The data for all energy points was combined together and the energy 
deposition of cosmic events was parametrized using the following function: 

fiE) = p ■ gi{E- Eo, a, 7?) + (1 - p) • gi{E; Eo + aa, /3a, 77?). (2) 
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Figure 7: Fit of the energy deposition for cosmic events with the function 
(H). Fit parameters: PI - number of events, P2 - p, P3 - i?o, P4 - a, P5 - 
r?, P6-a, P7-/3, P8-7 

The fit result is presented in Fig. |^. The obtained function was used as the 
energy deposition for background events: 

fcosmici-^) = fcosmici-^) — fcosmic{E), 

fco.rmc{E) = 0.69 • gi{E; 81.2, 9.0, -0.175) + (3) 

+0.31 • gi{E; 91.7, 22.8, -0.186). 

In order to check the stabihty of the energy deposition for cosmic events, 
all energy points were combined into 4 groups: 0.61-0.70 GeV, 0.71-0.778 
GeV, 0.780-0.810 GeV and 0.810-0.960 GeV. The energy deposition of the 
selected background events for each group was fitted with the function (|^). 
Results of the fit are demonstrated in Fig. ||. 

The probability to have no cluster in the calorimeter (zero energy depo- 
sition) for cosmic event tracks was also measured and found to be negligible 
— less than 0.1%. 
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(a) Group 0.61-0.70 GeV 
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(c) Group 0.780-0.810 GeV 
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(b) Group 0.71-0.778 GeV 
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(d) Group 0.810-0.960 GeV 



Figure 8: Fit of energy deposition of cosmic events with the function (g). 
Data is combined into 4 groups of energy points. Fit parameter: PI ~ number 
of events 
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(a) Fit of the positron energy depo- 
sition 



(b) Fit of the electron energy depo- 
sition 



Figure 9: Fit of the energy deposition of electrons and positrons at ccnter- 
of-mass energy 0.9 GeV. Fit parameters: PI - Pw, P2 - Eq, P3 - ct, P4 - ry, 
P5 - number of events, P6 - a, P7 - /3, P8 - 7 



2.2.6 Energy deposition of electrons and positrons 

The energy deposition for electrons and positrons was obtained from the data. 
Electrons and positrons can be separated from mesons by their relatively 
high energy deposition in the calorimeter. If we have one "clean" electron 
(positron) in the event (a particle with a proper momentum and high energy 
deposition), then we know that the second particle is a positron (electron), 
and therefore may be used for the energy deposition study. 

The test e+e~ e+e~ events were selected for the energy deposition 
study with the same selection criteria as for collinear events except for the 
following. 

1. The distance p from the vertex to the beam axis is less than 0.15 cm 
(instead of 0.3 cm). 

2. The average momentum of the particles is within the 10 MeV/c range 
around the beam energy. This requirement allows to decrease signifi- 



cantly a small admixture of e+e 



events. 



18 



Analysing the energy deposition in one cluster, the energy deposition in 
the other one is required to be between 0.92Eb — 100 MeV and 0.92Eb, where 
Eb is the beam energy. As for cosmic rays since the clusters are independent, 
the strict limit on one cluster energy does not affect the energy deposition 
distribution for the other cluster. 

For each energy point the energy depositions of the selected electrons and 
positrons were fitted with the following function: 

f{E) = (1-pw) ■ gi{E;Eo,a,'q) + ■ 9i{E;Eo + aa, /3cr,'^i]). 

After analysis of the fit results the following values of the parameters were 
found: 

a = -0.8, /3 = 3.0, 7=1.4. 

Again, as for cosmic rays, the probability to have no cluster was found to be 
negligible (less than 0.1%) for all energy points. Finally, the energy deposition 
of electrons and positrons was parametrized as: 

f+{E) EE f-^iE) EE U„{E), 

(4) 

fee{E) = {I -pn,)-gi{E-Eo, (7,7^) + p^ ■ gi{E; Eq - O.Sa, 3a, IAt)) , 

where p^, Eq, a and rj are free parameters of the fit. An example of the 
electron and positron energy depositions is shown in Fig. O. 



2.2.7 Energy deposition of muons 

In the centcr-of-mass energy range above 0.6 GeV, muons produced in the 
reaction e^e~ — *■ p,~^fJ-~ interact with the calorimeter as the minimum ion- 
ising particles, and therefore their energy deposition can be obtained from 
the detector simulation. The simulation of the e~^e^ — > i-J.^ p.^ events was 
performed for center-of-mass energies 0.6, 0.62, 0.64, 0.66, 0.68, 0.70, 0.74, 
0.75, 0.76, 0.77, 0.78, 0.79, 0.80, 0.82, 0.85, 0.90, 0.95, 1.0 GeV. 10000 events 
were simulated for each energy point. Simulated events were processed by 
the standard offline reconstruction program and the standard selection cri- 
teria for coUinear events were applied. Then the energy deposition of the 



simulated muons was corrected as described in sec. 2.2.4, 

The resulting energy deposition for each energy point was parametrized 
by the following function: 

fsiM^E) = Pi ■ [P2 ■ gi{E; Eo, ao, 1])+ 

+ (1 - P2) ■ 9i{E; Eo, /3o-o, 777)] + (1 - pi) ■ g{E; Ei,ai), 
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Figure 10: Simulated energy deposition of muons 
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where pi, p2, Eq, gq, Ei, gi, -q, j3 and 7 are the free fit parameters. 

In the simulation we have a perfectly calibrated calorimeter while the 
real data have an additional spread due to the uncertainty of the calibration 
coefficients. As a result, the real energy distribution is wider than that in the 
simulation. It could also be slightly shifted. To take that into account, the 
energy deposition distribution obtained from the simulation was modified to 
reflect the energy deposition of real muons: 

f+{E)^f;{E)^f^{E), 

UiE)^ f^upiE;k,a,), 

fMipiE;k,cr^) =pi • p2 ■ gi{E;Eo,&o,il) + 

+ {1-P2)- giiE;Eo,Pao,iv)] + (1 - Pi) ■ g{E; Ei, di), 
E = E/k, CTo = y/(TQ + 0-2, CTi = y/al +al. 

The values of pi, p2, Eg, txoi Ei, oi, rj, (3 and 7 are fixed according to 
the simulation, so this function has only two free parameters: k is the scale 
coefficient and <Tx is the additional energy spread. 

An example of the fit for one energy point as well as the energy dependence 
of Eq and ctq parameters are shown in Fig. pXl|. The k and ax parameters for 
all energy points are shown in Fig. fsk and |15|b respectively. 



2.2.8 Energy deposition of pions 

Pions produced in the reaction e^e~ — > tt+tt" may interact with the calorime- 
ter in two ways: as minimum ionising particles or through nuclear interaction. 

The energy deposition of the minimum ionising pions was obtained from 
the detector simulation in the same way as for muons. The simulation was 
performed for the center-of-mass energies 0.6, 0.62, 0.64, 0.66, 0.68, 0.70, 0.75, 
0.80, 0.85, 0.90, 0.95, 1.0 GeV. 10000 events were simulated for each energy 
point. Simulated events were processed by the standard offline reconstruction 
program and the standard selection criteria for coUinear events were applied. 
Then the energy deposition of the simulated pions was corrected as described 



in sec. 2.2A 



The resulting energy deposition for each energy point was parametrized 
by the same function as for muons: 

fsiiviiE) =pi ■ [p2 ■ giiE;Eo,ao,ri)+ 
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Figure 11: Simulated energy deposition of minimum ionising pions 
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+ (1 - P2) ■ 9i{E; Eo, f3<Jo,lv)] + (1 - Pi) ■ 9{E; Ei,ai), 



where pi, p2, Eq, aQ, cri, rj, [3 and 7 are the free fit parameters depending 
on the beam energy. This function describes weU the energy deposition only 
for those pions that do not stop in the calorimeter. 

As was done for muons, the modified simulated energy deposition was 
used for the parametrization of the energy deposition of real pions: 

fMipiE;k,cr^) =pi • p2 ■ gi{E;EQ,ao,T]) + 

+ (1 - P2) ■ 9i{E; Ea, /3ao, 777)] + (1 - pi) • g{E; Ei,<fi), 



E = E/k, CTo = + ^1 = V^f+af . 

The values of pi, p2, Eq, ao, Ei, ui, rj, (3 and 7 were fixed according to the 
simulation, so this function has only two free parameters: the scale coefficient 
k and the additional energy spread Ux- 

An example of the energy deposition of minimum ionising pions and the 
energy dependence of Eq and ctq parameters are shown in Fi g. The k and 
(Jx parameters for all energy points are shown in Fig. ^5|a and |l 5|b respectively. 

The energy deposition of the nuclear interacted pions cannot be simulated 
well. Therefore, in this case we have chosen to use an empirical parametriza- 
tion, which fits experimental data well: 

N 

g{E-0,cj)+Y.^g{E;^^i,a)-pl 



N ' 

where Eq is the beam energy, N is the number of Gaussian and a and pi.,,N 
are the free fit parameters. All Gaussian are normalized at the [0,-|-oo], 
therefore /^j is also normalized. This function appears as a set of equidistant 
Gaussians of the same width but different area. The parametrization with 
= 5 or 6 describes well the experimental data. 

The energy depositions of the minimum ionising tt"*" and 7r~ are the same 
while they are quite different for nuclear interacted tt"*" and Tr~ . Unlike elec- 
trons and muons, pions have a significant probability to have no cluster in 
the calorimeter and this probability is different for 7r+ and 7r~ . According to 
this, the overall pion energy deposition was parametrized with the following 
functions: 

f+iE)^p+-6iE) + il-p+)x (5) 
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X [pIiip ■ fliipiE; k, a^) + (1 - p+ ^p) • f;,i{E; a+, j3+ ,^)] , 
f-{E)^p, ■S{E) + {l-p^)x (6) 
X [pmip ■ fMip{E\k,a^) + (1 -p^/p) ' fNi{E\(r~,Pi,„j^)\ , 

where p^, P^npj a^, , p^ jy are free fit parameters. 

Pions produced in the (/)(i020) decay to Stt were used to test whether the 
functions and (^) describe the energy deposition of the real pions weh 
enough. The CMD-2 has collected about 20 pb~^ in the (yS- meson energy 
range. About 100000 "clean" 7r+ and 7r~ were selected from the completely 
reconstructed (p — s- tt^tt^tt^ events. The energy deposition of tt^ (''' ) with 
energies 345 < < 355 MeV and 425 < E^ < 435 MeV together with the 
fit are shown in Fig. |l2|a and [l2|3 (|l|b and |l^). 

2.2.9 Fit results 

The minimization of the likelihood function (^) was carried out for all energy 
points with the typical number of 10000 events per energy point with the 
following 27 free fit parameters: 

• global parameters 

ee I -''MM 

• Electron energy deposition parameters 

pw, Eq, (t, 77 

• Muon energy deposition parameters 

kMIP, 

• Pion energy deposition parameters 

TT + - + - 

^xi Po ' Po ' PmIP^ PMIPy 
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Figure 12: An example of the experimental energy deposition for pions of 
different energy. Pions were extracted from ip 7r"^7r~7r'^ data. The energy 
depositions of 7r+ and 7r~ were fitted with functions ^ and (ph respectively 
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Figure 13: Results of the minimization of the likelihood function (|l|). The 
energy dependences of the global parameters are shown 



The number of background events was fixed during the fit. The scale 
parameter for minimum ionising muons and pions was the same (listed about 
as Umip)- 

The results of the minimization are shown in Fig. p^-p^: the global pa- 
rameters are shown in Fig. [l^, the electron energy deposition parameters are 
shown in Fig.Jl4| and the muon and pion energy deposition parameters are 
shown in Fig. 15, 
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Figure 14: Results of the minimization of the likelihood function (|l|). The 
electron energy deposition parameters are shown 
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Figure 15: Results of the minimization of the hkehhood function (|l|). The 
muon and pion energy deposition parameters are shown 
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2.3 Form factor calculation 



X 



After the minimization of the hkehhood function, the pion form factor was 
calculated as: 

N 

\FJ^ = -211 X (7) 

Cr£ • (1 + See) (1 + "ee) £ee + O"^^ ' (1 + Sf,f,) (1 + a^^) 

cr^^ • (1 + Srr^) (1 + a^^) (1 - Ah) (1 - e^^ 

where A^7rTr/(A'ee + A'^^) was obtained as a result of minimization, cr^ are the 
Born cross-sections, Sa are the radiative corrections, Sa are the efficiencies, 
including trigger and reconstruction efficiencies, aa are the corrections for the 
experimental resolution of the O angle measurement. Ah is the correction for 
pion losses due to nuclear interactions, A/^ is the correction for pion losses 
due to decays in flight and Aa^r is the correction for the misidentification of 
7r"'"7r~7r'^ events as a 7r"''7r~ pair. 

Together with the pion form factor, the cross-section of e+e^ annihilation 
to TT+TT^ was calculated as 



and the luminosity integral as 
j Ldt= - 



Ami 



3/2 



• \FJ\ 



Nee + N^^ 



(1 + 5ee){l + aee)£ee + O^n ' (1 + 5^,^,){l 



2.3.1 Radiative corrections 

The radiative corrections are shown in Fig. ^ The calculation of the radia- 
tive corrections for e^e'^ e~^e~ events was based on The estimated 
systematic error is ~ 1%. Radiative corrections for e+e^ tt+tt" 
were calculated according to and The estimated systematic error 

is 0.2 — 0.5%. The contribution from the lepton and hadron vacuum po- 
larization is included in the radiative corrections for e~^e~ e'^ e~ , ^'^ fj,~ , 
but excluded from the radiative correction for e+e^ — > tt+tt^. Some details 
about the radiative correction calculation could be found in |^ . 

The radiative correction for e'^e~ — > tt'^tt~ depends on the energy be- 
haviour of the e+e^ — > tt+tt^ cross-section itself. In order to take that into 
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Figure 16: Radiative corrections for e+e — > e+e ,7r+7r events 



account, the calculation of this correction was done via a few iterations. At 
the first iteration, the existing |f^(s)|^ data were used for the calculation of 
the radiative correction. Using the calculated correction, we obtain |-F7r(s)| • 
At subsequent iterations, ji^Trls)!^ data obtained from the previous step, were 
used. It was found that after 3 iterations |-F7r(s)| values become stable. 

2.3.2 Trigger efficiency 

The trackfinder (TF) signal, neutral trigger (NT) signal and the calorimeter 
"OR" (CSI) trigger signal have been used for detector triggering during 1994- 
1995 runs. The trackfinder uses only the drift chamber and the Z-chamber 
data for a fast decision whether there is at least one track in the event. 
The neutral trigger analyses the geometry and energy deposition of the fired 
calorimeter rows. The calorimeter "OR" signal appears when there is at least 
one triggered calorimeter row. 

Two different trigger settings have been used during 1994-1995 runs. 
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For 810(2)-960 MeV energy points the trigger was "(TF.and.CSI).or.NT" 
while for other energy points the trigger was "TF.or.NT". For analysis we 
use only those events where the trackfinder has found the track. If the track- 
finder efficiency is etf and the CSI signal efficiency is ecsi, the overall 
trigger efficiency is etf ■ £csi for 810(2)-960 MeV energy points and etf for 
610-810(1) MeV energy points. 

For the measurement of the trackfinder efficiency e+e^ — > e+e^ events 
have been selected using only the calorimeter data (see the next section for de- 
tails). Additionally, the positive decision of the neutral trigger was required. 
Analysing the probability to have also a positive decision of the trackfinder, 
its efficiency for e^e^ e+e^ events can be calculated. The result is shown 



in Fig. 17(a). All types of coUincar events are similar from the trackfinder 
point of view — they have two coUinear tracks and the track momenta are 
very close for different types of events. It is seen from the picture that the 
TF efficiency is high enough and that there is no visible energy dependence 
of the efficiency. Therefore we assume that the trackfinder efficiency is the 
same for all types of coUinear events and hence cancels in (|^). 

For the measurement of the CSI signal efficiency the fact that there are 
two clusters in the coUinear event have been used: the CSI efficiency for a 
single cluster has been measured and was used for the calculation of the CSI 
efficiency for the event. The coUinear events have been selected using the 
standard cuts. Then, the probability to trigger the CSI signal was measured 
for one cluster as a function of its energy. Using the measured probability 
and our knowledge about the energy deposition of particles of different types, 
the efficiency to trigger the CSI signal was calculated for different ty pes of 
events. Results for 810(2)-960 MeV energy points are shown in Fig. 17(b) . 
For simplicity, the values of 99.5% for pious and 99.2% for muons for the CSI 
efficiency were used for all 810-960 MeV energy points. 



2.3.3 Reconstruction efficiency 

The reconstruction efficiency was measured for e+e~ —^ e~^e~ events using 
the experimental data: events were selected using only the calorimeter data, 
and then the probability to reconstruct two coUinear tracks was measured. 
Since the selection criteria for coUinear events are based on tracking data 
only, the measured probability is the reconstruction efficiency. 

The following selection criteria based only on the calorimeter data, were 
used. 

1. There are exactly two clusters in the calorimeter. 
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Figure 17: Trigger efficiency for collinear events 
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Figure 18: Reconstruction efficiency for all energy points 



2. There is a hit in the Z-chambcr near each cluster. This requirement 
selects the clusters produced by the charged particle. 

3. The energy deposition of both clusters is between (0.82 • Eb — 40) and 
(0.82 • Eb + 50) MeV. This is typical for e+ and e", but not for other 
particles. 

4. The clusters are collinear if one takes into account the particle motion 
in the detector magnetic field: 



|7r-(ei + e2)| <o.i, 

|<^i-<^2||-2arcsin(^-3^ 



sin 



< 0.1, 



where O and ip are the polar and azimuthal angles of the cluster, R 
is the effective calorimeter radius (45 cm) and B is the magnetic field 
(typically 10 kGs). 

5. The event was triggered by the trackfinder. Other triggers may be 

present. 

The selected sample contain only e+e~ e+e~ events with negligible (< 
0.5%) background. If the number of selected events is Nq and the number of 
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events, which in addition satisfy the standard selection criteria for colhnear 
events is A^i, then the reconstruction efficiency is 



The measured reconstruction efficiency is shown in Fig. ^ (filled dots). The 
efficiency drops down for beam energies below 400 MeV. This additional 
inefficiency is completely explained by the lower precision of the angle 
measurement for these energy points. The empty dots in the same plot rep- 
resent the reconstruction efficiency for the same selection criteria for collinear 
events except for AG. 

As in studying the trigger efficiency, all types of collinear events are very 
similar from the point of view of the tracking reconstruction code. Therefore 
we assume that the reconstruction efficiency is the same for all types of 
collinear events and hence it cancels in (^. 

2.3.4 Correction for the limited accuracy of the polar angle mea- 
surement 

There is a special type of correction that arises since there is finite exper- 
imental resolution in the measurement of the Q angles of tracks. Since we 
require the measured, rather than real, Qavr = [©i -I- (tt — 62)]/2 to be 
within [Qmin, TT — ©mm] interval, the visible cross-section is slightly different 
from 

— • 27rsinede. 

0^.^ ail 

The corresponding correction, denoted a in (0), was calculated as 

1 + J de J de' ■ fie')p{&,e) / J de-/(e), 

e,„,„ ' e„i„ 

where /(6) is the theoretical 6-distribution and p(6',9) is the detector 
response function taken as 



p(e',e)--=^cxp 



(9' - 6)' 
2a| 



where is the experimental resolution of Qavr measurement. The resolution 
and the corresponding corrections for different types of collinear events for 
all energy points are shown in Fig. 
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Figure 19: The experimental resolution of the average polar angle measure- 
ment and the corresponding correction to cross-sections 



2.3.5 Correction for pion losses due to nuclear interactions 

Unlike electrons and muons, a small fraction of pions has nuclear interactions 
in the beam pipe and therefore is lost from the selected samples of collinear 
events. The corresponding correction was calculated under the assumption 
that the event is lost from the collinear event sample if at least one pion had an 
inelastic nuclear interaction in the beam pipe or drift chamber inner support. 
The calculations were done separately for cross-sections obtained from two 
packages available in GEANTjlJ] for simulation of nuclear interactions — 
and GHEISHA||l6|]. The resulting correction is shown in Fig. 
has been shown that FLUKA is more successful in describing 



FLUKAflS 
1^. It 

experimental data about nuclear interactions of pions in our energy range, 
and therefore the results obtained from FLUKA were used for the cross- 
section correction while the difference between two calculations was added to 
the systematic error. 



2.3.6 Correction for pion losses due to decays in flight 

A small fraction of pions decays in flight inside the drift chamber and the 
corresponding event may not be recognized as collinear if the decay angle 
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Figure 20: The probability to lose the e+e~ — > tt+tt" event from the 
sample of collincar events due to nuclear interactions of pions with the beam 
pipe or drift chamber inner support 
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Figure 21: The probability Ad to lose the e+e — > tt+tt event from the 
sample of coUinear events due to pion decay in flight 
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Figure 22: The A^^^ contribution to the pion form factor from a; — )■ Stt 



between the pion and secondary muon is large enough. Depending on the 
beam energy 2-4% of the e~^e~ tt+tt" events have a decayed pion. But 
the maximum decay angle for the analysed energy range is small — 80-150 
mrad depending on the pion energy — and therefore more than 90% of such 
events are reconstructed as a pion pair. The probability to lose a pion pair 
due to the pion decay in flight was found from simulation and is shown in 
Fig. |1[ 

2.3.7 Correction for the background from uj ^ 3tt 

There is a small probability to identify the tt+tt^tt" final state as a tt+tt^ pair. 
Therefore the measured pion form factor is slightly bigger than it should be. 
This effect is significant only in the narrow energy range near the w-meson. 
The corresponding probability was estimated from simulation and was found 
to be 5-10~^ for Qmin = 1-0 and 4-10"'^ for Qmin — 1-1 rad. The contribution 
to the pion form factor A^t^ was calculated assuming the w-meson parameters 
from and is presented in Fig. 
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2.4 Systematic errors 

The overall systematic error is estimated to be 1.5% for 0.780 and 0.784 GeV, 
1.7% for 0.782, 0.84-0.87 and 0.94 GeV and 1.4% for other energy points. The 
main sources of systematic error are summarized in Table |^ and discussed in 
detail below. 



Source 


Estimated value 


Events separation 


0.6% 


Energy calibration of collider 


0.1% (0.5% for 390,392; 1% for 391) 


Fiducial volume 


0.5% 


Trigger efficiency 


0.2% (1% for 420-435, 470) 


Reconstruction efficiency 


0.3% 


Hadronic interactions of pions 


0.4% 


Pions decays in flight 


0.1% 


Radiative corrections 


1% 


Total 


1.4% 




(1.5% for 390,392) 




(1.7% for 391,420-435,470) 



Table 2: Main sources of systematic errors 



Event separation. Special studies were performed to estimate the system- 
atic error due to event separation. 

1. The separation was done by minimization of two independent likelihood 
functions. The "standard" likelihood function was described in detail 
in the previous sections. The second likelihood function was different 
in the following ways: different separation of background events; differ- 
ent energy deposition functions for minimum ionising muons and pions 
and for the nuclear interactions of pions (a single Gaussian was used 
instead of the set of Gaussians). The relative difference between two 
minimizations results is presented in Fig. The error bars represent 
the statistical error. It is seen that there is no systematic shift and the 
difference is well below the statistical error. 

2. The complete detector simulation of e+e~ e+e~(7), tt+tt" 
events in the GEANT environment was done for 10 energy points: 
0.61, 0.62, 0.64, 0.66, 0.70, 0.778, 0.782, 0.786, 0.860 and 0.960 GeV. 
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(b) The relative difference between simulated and recon- 
structed form factor values for the complete detector sim- 
ulation of coUineaj: events 



Figure 23: Estimation of the systematic error of events separation procedure. 
The relative difference A|F^p/|F^p is shown with the error bars representing 
statistical error 
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100000 of e^e~ —* 6^6^(7) events and the corresponding number of 
e+e^ TT+TT^ events were simulated for each energy point. 

The event separation was done by minimization of the same hkehhood 
function as for real events. The relative difference between simulated 
and reconstructed form factor values is shown in Fig. ^3|b. There is no 
systematic shift for all energy points except for the first three (0.61, 
0.62 and 0.63 GeV). 

3. For some energy points there is a small difference between energy de- 
positions of electrons and positrons, which comes mainly from a dif- 
ference in the energy calibration of different parts of the calorimeter. 
The corresponding systematic error was estimated from the simula- 
tion. 100000 of e+e" e+e~ events and the corresponding number of 
e^e~ — > TT^TT" events were simulated for 7 energy points. The 

energy deposition of particles was generated according to the energy 
deposition functions described in the previous sections, but with differ- 
ent parameters for e+ and e~ . Then the event separation was done by 
minimization of the same likelihood function as for the real events. The 
systematic shift of about 0.5 — 0.8% was observed between simulated 
and reconstructed form factor values. 

We estimate the average systematic error of event separation to be less 
than 0.6%. For the first three energy points the systematic error is larger 
(4%, 2% and 1% for 0.61, 0.62 and 0.63 GeV energy points respectively), but 
since the statistical error is a few times larger we neglect this fact. 



Energy calibration of the collider. For almost all energy points the 
beam energy was measured using the resonance depolarization technique |l9|] . 
The systematic error in beam energy measurement does not exceed 50 keV. 
The corresponding systematic error in the pion form factor does not exceed 
0.1% for all energy points, except for the narrow p ~ to interference energy 
range, where it is estimated to be less than 0.5% for 390 and 392 MeV energy 
points and less than 1% for 391 MeV energy point. 



Fiducial volume (accuracy of Qmin measurement). The Z coordinate 
of the track intersection with the Z-chamber is measured with better than 1 
mm systematic accuracy. That allows, after the calibration of the drift cham- 



ber by the data from the Z-chamber 1 20 , to have systematic accuracy in polar 



angle measurement better than 2.5 mrad. The corresponding systematic er- 
ror in the pion form factor is 0.5%. Since the pion form factor was measured 
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independently for two Qmin values , we can estimate this systematic error 
by the average form factor deviation (l — |e„i„=i.o /I^ttP |e„i„=i.i )• 
This deviation is shown in Fig. ^ It is equal on the average to — 0.1%±0.3% 
and is consistent with the expected statistical fluctuation. 



U. 2 - 




800 900 

Center-of-mass Energy, MeV 



Figure 24: The relative difference in the form factor values, measured for 
©mm — 1-0 and Qmin = 1-1 rad. The errors correspond to statistical fluc- 
tuation of the number of events in 1.1 < < 1.0 region. In other words, 
the errors show the expected statistical difference for the form factor values 
measured for Qmin = 1.0 and ©mm = 1.1 



Trigger efHciency. For most energy points the trackfinder efficiency is 
high (> 98%), and as it was discussed above, all types of coUinear events are 
similar from the trackfinder point of view. Therefore it is assumed that the 
corresponding systematic error is negligible. This statement is questionable 
for 420-435 and 470 MeV energy points where the TF efficiency is lower 
than usual. The Csl efficiency is high for all energy points (> 99%) and 
was measured for all types of collinear events. Therefore the corresponding 
systematic error is negligible. The combined systematic error is estimated to 
be less than 0.2% for all energy points except 420-435 and 470 MeV, where 
it is estimated to be less than 1%. 

Reconstruction efficiency. The reconstruction efficiency was measured 
only for Bhabha events, but it was found that it is uniform over ip and 
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and that it comes mostly from errors in the track recognition algorithm. 
Therefore it is reasonable to assume that the reconstruction efficiency is the 
same for all types of collinear events. The efficiency drop for energy points 
below 800 MeV caused by the lower precision of 8 angle measurements is 
also the same for all types of collinear events. The corresponding systematic 
error is estimated to be less than 0.3%. 

Correction for the nuclear interaction of pions. The estimated ac- 
curacy of the correction calculation is about 20%, determined by the cross- 
sections precision. The estimated contribution to this correction from nuclear 
interaction of pions inside the drift chamber and from elastic pion scattering 
by a large angle is about 30% of the correction. Therefore the corresponding 
systematic error is less than 0.4%. 

Correction for the pions decays in flight. The systematic error of this 
correction is estimated to be less than 0.1%. 

Radiative corrections. The systematic error in radiative corrections is 
dominated by the error in radiative corrections for e+e^ — > e^e~, where it 
is estimated to be less than 1%. 

3 A fit of the pion form factor 

The pion form factor data, obtained from the described analysis, is summa- 
rized in Table ^j. Only statistical errors are shown. 

ft is well known that to describe the pion form factor the higher resonances 
p(1450) and /3(1700) should be taken into account in addition to leading 
contribution from p(770) and w(782) H, |2l[. At the same time it was shown 
I p^ that the experimental data below 1 GeV is well described by the model 
based on the hidden local symmetry, which predicts a point-like coupling 
77r"^7r~. Here we use both approaches to fit the data. Since we fit the pion 
form factor in the relatively narrow energy region 0.61-0.96 GeV, only one 
higher resonance p(1450) is taken into account. 

3.1 The Gounaris-Sakurai (GS) parametrization 

This model is based on the Gounaris-Sakurai parametrization of the p-reso- 
nance. Its modifications were used for parametrization of all previous e^e~ 
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2E (McV) 




CTttti- (nb) 


— iT 


610.50 


426 ± 50 


333.6 ±38.9 


8.14 ±0.95 


620.50 


829 ± 60 


378.0 ±26.8 


9.41 ±0.67 


630.50 


948 ± 64 


422.8 ±27.8 


10.74 ±0.71 


640.51 


1005 ± 63 


436.3 ±26.6 


11.31 ±0.69 


650.49 


851 ±55 


492.2 ±30.2 


13.02 ±0.80 


660.50 


1288 ±65 


549.4 ±26.2 


14.83 ±0.71 


670.50 


1513 ±74 


602.6 ±27.8 


16.60 ±0.77 


680.59 


1654 ±77 


697.9 ±30.0 


19.63 ±0.84 


690.43 


2257 ± 87 


717.8 ±25.3 


20.59 ±0.73 


TOO. •12 


;]411 ± 108 


84.1.1 ± 2;]. 9 


24.7()± 0.70 


710.47 


4196 ± 131 


946.0 ± 26.3 


28.28 ± 0.79 


720.25 


3340 ± 113 


1018.3 ±30.2 


31.06 ±0.92 


730.24 


2764 ± 108 


1111.9 ±37.7 


34.61 ± 1.17 


740.20 


3138 ± 116 


1181.5 ±37.7 


37.53 ± 1.20 


750.28 


4290 ± 143 


1315.9 ±37.3 


42.66 ± 1.21 


760.18 


3745 ± 136 


1337.8 ±41.0 


44.24 ± 1.36 


764.17 


6241 ± 168 


1301.8 ±29.2 


43.40 ±0.97 


770.11 


4805 ± 151 


1301.1 ±34.3 


43.90 ± 1.16 


774.38 


3729 ± 127 


1250.5 ±35.8 


42.55 ± 1.22 


778.17 


3897 ± 135 


1313.6 ±38.0 


45.04 ± 1.30 


780.17 


3993 ± 136 


1253.9 ±35.7 


43.17 ± 1.23 


782.23 


9709 ± 204 


1080.7 ± 19.2 


37.36 ±0.66 


784.24 


4946 ± 143 


986.9 ±24.3 


34.25 ±0.84 


786.04 


2709 ± 108 


869.6 ±30.4 


30.29 ± 1.06 


790.10 


4948 ± 149 


886.3 ±23.3 


31.12 ±0.82 


794.14 


4(iG(i ± 14.5 


83G.-5± 2;].() 


29.61 ± 0.81 


800.02 


4701 ± 139 


819.6 ±21.2 


29.35 ±0.76 


810.14 


6042 ± 152 


699.2 ± 15.5 


25.55 ±0.57 


809.98 


6648 ± 163 


696.5 ± 15.0 


25.44 ±0.55 


820.02 


2606 ± 99 


648.2 ±21.8 


24.15 ±0.81 


829.97 


2209 ± 90 


542.7 ± 19.9 


20.61 ±0.76 


839.10 


1687 ±77 


450.3 ± 18.6 


17.41 ± 0.72 


849.24 


1194 ±51 


362.2 ± 14.5 


14.28 ±0.57 


859.60 


1288 ± 63 


364.3 ± 16.6 


14.65 ±0.67 


continued on the next page 



Table 3: The experimental data from the CMD-2 detector 
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2E (MeV) 




(Jtttt (nb) 




869.50 


1581 ±71 


274.8 ± 11.6 


11.26 ±0.48 


879.84 


544 ± 42 


260.6 ± 19.1 


10.89 ±0.80 


889.72 


1610 ± 70 


203.4 ±8.4 


8.66 ±0.36 


900.04 


1885 ± 76 


182.0 ±7.0 


7.90 ±0.30 


910.02 


1352 ±65 


156.4 ±7.3 


6.92 ±0.32 


919.56 


1149 ±60 


135.6 ±6.9 


6.11 ±0.31 


930.11 


748 ± 48 


127.8 ±8.0 


5.87 ±0.37 


942.19 


1272 ±61 


116.1 ±5.4 


5.45 ±0.26 


951.84 


1220 ±66 


96.3 ±5.1 


4.60 ±0.24 


961.52 


1097 ±60 


91.6 ±4.9 


4.45 ±0.24 



Table 3: The experimental data from the CMD-2 detector 



data [|j and r decay data To describe the data, the p{770), /9(1450) 

contributions and p — lo interference were taken into account: 



BW 



GS 

p(770) 



is) 



l + (5BW^(.s) 
1 + 5 



+ /3BW«fi45o)(s) 



1 + 13 



(8) 



The p(770) and p(1450) contributions are taken according to the Gouna- 
ris-Sakurai model lESl : 



GS 



Ml{\ + d-Tp/Mp) 



(9) 



where 



Ml 



dh 
ds 



dh 



pUs) {h{s) - hiMD) + {Ml - s)pI{MI) — 



s=Ml 



2m^ 



= HMl 



,2(M2) 2M2 



27rM2' 
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(10) 
(11) 

(12) 
(13) 



The energy dependent width is 





3 


Ml 






s 



1 1/2 



(14) 



The normahzation BW^(^2 •) (0) = 1 fixes the parameter d: 



.2 Mp + 2p,(M2) 



In- 



mlMp 



2m^ 



(15) 



The simple Breit-Wigner parametrization ( p5[) is used for the ^(782) con- 
tribution. 

In order to extract T{p — > e^e~) we relate the pion form factor at the p 
mass to the VMD form factor: 



From this one can obtain 



M2 - s - iMpTp 



non — resonant 
contribution 



Ml{l + d-Tp/Mp) 



Using well known VMD relations 



iira^ 2 



3M' 



V 



V — ^7T+7r 



and assuming that Tp^^+^- = Tp one obtains 



9MpFp (1 + ^)2(1 + /3)2 



(16) 

(17) 

(18) 
(19) 

(20) 



A similar approach is used for the calculation of Br{uj — > tt+tt ). One 
can relate the pion form factor at the uj mass to the VMD form factor: 



=M2 



Ml-s- iM^T^ 



non — resonant 
contribution 



(21) 
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From this one determines g^. 



<5.MMBW°f„o)(M^y 



'yyuJTTTT 



(22) 



Using (|18| ) and one derives 



Br{bJ TT+TT ) 



9M^r 



BW 



p(770)(M^ 



(1+ (5)2(1 +/3)2 



(23) 



3.2 The Hidden Local Symmetry (HLS) parametriza- 
tion 

In the Hidden Local Symmetry (HLS) model j22[ the p-meson appears 

as a dynamical gauge boson of a hidden local symmetry in the non-linear 
chiral Lagrangian. This model introduces a real parameter a related to the 
non-resonant coupling 77r+7r^. The resulting pion form factor for the HLS 
model is 

F.(^)--i + l + i-BW,(.).i±i5^, (24) 



where 



BWy(s) 



iMvTv{s)' 



(25) 



The energy dependent width Tp{s) is the same as for the GS parametriza- 
tion and is given by (p^. 

Using the similar approach as for the GS parametrization, we obtain 



r 

^ p — >e 



Br{u! — > TT^TT ) = 



2a'pl (Mg) 
OMpLp 



1 2 



2(1 + 5) 



2BWp(M2) 



(1 + <S)2 



(26) 
(27) 



3.3 Result of fit 

The fit to data was done by the minimization of the following function: 



X 



E 

i=1...44 



3 = 1. ..Np 



{Pj-POjY 



(28) 
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GS model 


HLS model 


Mp, MeV 


775.28 ±0.61 ±0.20 


774.57 ±0.60 ±0.20 


Tp, MeV 


147.70 ± 1.29 ±0.40 


147.65 ± 1.38 ± 0.20 


Br{uj TT+TT-), % 


1.31 ±0.23 ±0.02 


1.32 ±0.23 ±0.02 


r(p-> e+e-), keV 


6.93 ±0.11 ±0.10 


6.89 ±0.12 ±0.10 


/3 (GS) 


-0.0849 ±0.0053 ±0.0050 




a (HLS) 




2.381 ±0.016 ±0.016 




0.77 


0.78 



Table 4: The results of fit of the CMD-2 (94,95) pion form factor data by GS 
and HLS models 

where I^VlLpl'^i) ^^'^ \^Tr\theori^i) experimental and theoretical val- 

ues of the pion form factor at the i-th energy point, is the experimental 
error at the i-th energy point, Np is the number of model parameters, Pj is 
the j-th model parameter, poj is the expected value of j-th model parameter 
and (Tj is the estimated error of j-th model parameter. The model parameters 
are the following (all values are taken from PDG'98 p8[). 

1. The uj meson mass: poi = 781.94 MeV, ai = 0.12 MeV. 

2. The uj meson width: po2 = 8.41 MeV, (72 = 0.09 MeV. 

3. The u! meson leptonic width: po3 — 0.60 keV, 1T3 — 0.02 keV. 

4. The /9(1450) meson mass: po4 = 1465 MeV, CT4 = 25 MeV. 

5. The p(1450) meson width: po5 = 310 MeV, CT5 = 60 MeV. 

For the GS fit the number of model parameters is Np ~ 5, while for the HLS 
fit Np = 3. 

The parameter /3 in (^ was assumed to be real. To determine a phase of 
p — UJ mixing, the parameter 5 in (||) and (^J) was assumed to be complex 
during the fit. However, it turns out that the phase was always consistent 
with zero. Therefore 5 was assumed to be real for the final fit. 

The results of the fit of the CMD-2 experimental data are shown in Fig. ^ 
and are summarized in Table ^ The first error is statistical, the second error 
is systematic. In order to estimate a systematic error, the fit was repeated 
with all data points shifted up and down by one systematic error. 
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Figure 25: Fit of the CMD-2 (94, 95) pion form factor data according to 
GS and HLS models. Both theoretical curves are presented, but they are 
indistinguishable 
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3.4 Hadronic Contribution to (g-2) 



We'll now estimate the implication of our results for the value and uncer- 
tainty of a^"-'^ - the hadronic contribution to (g-2)^. To this end we'll choose 
the cm. energy range near the p-meson peak (from 630 to 810 MeV) which 
gives the dominant contribution to the muon anomaly and where good data 
are available from previous measurements by the OLYA and CMD groups Q 
as well as by the DM1 detector Table ^ presents results of our calcula- 
tions of ajj"'' performed by direct integration of the experimental data over 
the energy range above. The calculation used the following procedure: first 
we integrated the data of each group over the chosen energy range. Since 
measurements at different energy points are independent, the squared sta- 
tistical error of the integral can be obtained by summing statistical errors 
squared of each energy point. After that the systematic uncertainty which 
is believed to be an overall normalization uncertainty is added quadratically. 
Contributions of separate groups are subject to weighted averaging taking 
into account if necessary a scale factor [ |l8| . Such a procedure implies an 
assumption that systematical uncertainties are uncorrelated since they refer 
to different measurements. 



Data 




Total error, lO^^'' 


Old 


284.2 ± 3.8 ± 7.1 


8.1 


New 


292.5 ± 2.2 ± 4.1 


4.7 


Old-hNew 


290.3 ± 1.9 ± 3.6 


4.1 



Table 5: Hadronic contributions to (g-2)^, coming from the cm. energy range 
near the p-meson peak. 

The first line gives the average of three independent estimates based on 
the data of OLYA, CMD and MDl while the second one presents the cor- 
responding value for the CMD-2 data. The third line of the Table is the 
average based on the four independent estimates above. For convenience, 
we list separately the statistical and systematic uncertainties in the second 
column while the third one gives the total error obtained by adding them 
quadratically. One can see that the estimate based on the CMD-2 data is in 
good agreement with that coming from the old data. Note also the signifi- 
cant improvement of both statistical and systematic uncertainties compared 
to the previous measurements. The combined result (old plus new data) is a 
factor of two more precise than before. 
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4 Conclusion 



A new measurement of the pion form factor in the center-of-mass energy 
range 0.61 — 0.96 GeV is presented. The statistical error is about the same 
as for aU previous e'^e~ data, but the systematic error is a few times smaller. 
The results are based only on part of the experimental data taken by CMD- 
2. There are several reasons why this data set was treated as independent. 
First, a different approach in data analysis is required for the data taken 
in the center-of-mass energy range below 0.6 GeV and above 1.0 GeV (runs 
4 and 5 in table |l]). Second, the data taking conditions were significantly 
different in 1998 (run 6 in table [|), when the data was taken in the same 
0.61 — 0.96 GeV energy range. As a result, the systematic errors for these 
runs will be different than for the data presented here. Therefore we prefer 
to present separate results for different groups of runs. 

The fit of the pion form factor based on the Gounaris-Sakurai and the 
Hidden Local Symmetry parametrizations was performed. For the final values 
we prefer the Gounaris-Sakurai parametrization as it was traditionally used 
for the determination of the p-meson parameters (see ||] and [^). The 
difference between the results of the fit in two models is always smaller than 
the experimental error (Table |^), therefore the "model" error is not specified. 
We give as the final results the following, where the first error is statistical 
and the second one is systematic: 



Mp(MeV) = 775.28 ±0.61 ±0.20, 
rp(MeV) 147.70 ± 1.29 ±0.40, 

r(p-^e+e-)(keV) = 6.93 ± 0.11 ± 0.10, 

Br{uj ^ n+TT-) = (1.31 ±0.23)%. 



(29) 



It has also been shown that the improvement of the experimental preci- 
sion, particularly of the systematic uncertainties, can be crucial for the high 
precision calculation of the hadronic contribution to the muon anomaly. 

With the progress of analysis of the rest available data we hope to reduce 
the systematic error for the data presented here by additional factor of two. 
Important part of this improvement is the development of the new approach 
to the radiative corrections calculation with the ultimate goal to reduce the 
corresponding systematic error to (0.3 — 0.5)% level. Using the high statistics 
taken in 1998, it will be possible to perform the detailed analysis of small 
systematic effects. 
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